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Abstract. We develop a scheme for the extraction of the properties of the scalar mesons /o(600), /o(980), 
and ao(980) from lattice QCD data. This scheme is based on a two-channel chiral unitary approach with 
fully relativistic propagators in a finite volume. In order to discuss the feasibility of finding the mass and 
width of the scalar resonances, we analyze synthetic lattice data with a fixed error assigned, and show 
that the framework can be indeed used for an accurate determination of resonance pole positions in the 
multi-channel scattering. 
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1 Introduction 

The application of lattice QCD techniques to obtain ha- 
dronic spectra and resonances is catching up and looks 
very promising for the near future [HHIMlIMlElillilini 
[TT | |T ^ 3, 14, 15. lSl[T71|T51[Tgi[^ . However, as it is well known 
resonances do not correspond to isolated energy levels in 
the (discrete) spectrum of the QCD Hamiltonian mea- 
sured on the lattice, and an additional effort is needed 
to extract the parameters of a resonance (its mass and 
width) from the "raw" lattice data. In case of clastic scat- 
tering, the pertinent procedure is well known under the 
name of Liischer framework [^Tll22| . In this framework, 
for a system described by a given quantum-mechanical 
Hamiltonian, one relates the measured discrete value of 
the energy in a finite volume to the scattering phase shift 
at the same energy, for the same system in the infinite 
volume. Consequently, studying the volume-dependence 
of the discrete spectrum of the lattice QCD gives the en- 
ergy dependence of the elastic scattering phase shift and 
eventually enables one to locate the resonance pole posi- 
tions. 

Recently, the Liischer approach has been generalized | 
to the case of multi-channel scattering. This was done in 
Refs. [23] on the basis of potential scattering theory, while 
the Refs. [^lU^ use non-rclativistic effective field theory 
(EFT) for this purpose. As discussed in Ref. [5S], there is 
a fundamental difference between Liischer approaches in 
the elastic and the inelastic cases: whereas in the former, 



one aims at the extraction of a single quantity (the scat- 
tering phase shift) from the single measurement of the en- 
ergy level, there will be, e.g., three different observables at 
a single energy in case of two-channel scattering (conven- 
tionally, the two-channel S'-matrix is parameterized by two 
scattering phases and the inelasticity parameter) . In order 
to circumvent the above problem, in Ref. [15] it has been 
proposed to impose twisted boundary conditions (b.c). 
Performing the measurements of the spectrum at differ- 
ent values of the twisting angle provides one with data, 
which are in principle sufhcicnt to determine all S'-matrix 
elements in the multi-channel scattering independently. 
Alternative proposals, e.g., using asymmetric lattices, will 
be briefly considered in the present paper. 

In this paper we consider a framework for the extrac- 
tion of the scalar resonance parameters, based on unita- 
rized Chiral Perturbation Theory (UCHPT). In the in- 
finite volume, this model is very successful, and repro- 
duces well the TTTr/inj and KK data up to 1200 MeV. 
The resonances /o(600), /o(980), ao(980) are generated dy- 
namically from the coupled channel interaction [2M77[f251 
^ . In this paper, we consider the same model in a finite 
volume to produce the volume-dependent discrete energy 
spectrum. Reversing the argument, one may fit the pa- 
rameters of the chiral potential to the measured energy 
spectrum on the lattice and, at the next step, determine 
the resonance locations by solving the scattering equations 
in the infinite volume. 
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In this work wc address two main issues. The first one 
is the use of fully relativistic propagators in the effective 
field theory framework in a finite volume. It is demon- 
strated that this allows one to avoid the sub-threshold sin- 
gularities which are inherent to the Liischer approach. In 
case of the multi-channel scattering, these sub-threshold 
singularities may come dangerously close to the physical 
region or even sneak into it. We also study the numerical 
effect of the use of relativistic propagators on the finite- 
volume spectra and on the extraction of the physical ob- 
servables from the lattice data. 

Second, we discuss in detail the analysis of "raw" lat- 
tice data for the multi-channel scattering. To this end, we 
supplement lattice data by a piece of the well-established 
prior phcnomenological knowledge that stems from UCHPT 
in order to facilitate the extraction of the resonance pa- 
rameters. In particular, it will be shown that, with such 
prior input, e.g., the extraction of the pole position from 
the data corresponding only to the periodic b.c, is indeed 
possible. 

In the present paper, in order to verify the above state- 
ments, we shall analyze "synthetic" lattice data. To this 
end, we produce energy levels by using UCHPT in a fi- 
nite volume, assume Gaussian errors for each data point, 
and then consider these as the lattice data, forgetting how 
they were produced (e.g., forgetting the parameters of the 
effective chiral potential and the value of the cutoff). In 
the analysis of such synthetic data, we shall test our ap- 
proach, trying to establish resonance masses and widths 
from the fit to the data. Note that our procedure is sim- 
ilar to the one used previously in Ref. [3D] for the case 
of the Z\-resonance in -kN elastic scattering. Recently, the 
problem of extracting continuum quantities from synthetic 
lattice data has also been also discussed in Ref. [SI] in 
the context of meson-meson interaction with charm and 
the Z)s*o(2317) resonance. The extraction of continuum 
quantities, at higher energies, may also be performed us- 
ing dynamical coupled-channels approaches as shown in 
Ref. [SI]. 

The paper is organized as follows. In Sec. [5] we develop 
the relativistic effective field theory formalism in a finite 
volume. In section [3] we elaborate our method in coupled 
channels and discuss the use of twisted b.c. and asymmet- 
ric lattices for the extraction of the multi-channel reso- 
nances. Finally, in section H] we give a detailed discussion 
of the fit by using UCHPT language in a finite volume. 
Section [^ contains our conclusions. 



2 Formalism 

2.1 General setting 

Below, we briefly present our approach to the finite vol- 
ume problem within the setting of UCHPT. Unless stated 
otherwise, we shall restrict ourselves to the case with total 
isospin 7 = (the case with 1=1 can be considered anal- 
ogously and is briefly discussed in section [53]). We thus 
start by using the chiral unitary approach in two channels, 
TTTT and KK, with the potential obtained from the lowest 



order chiral Lagrangians [33] ■ We follow the approach of 
Ref. [ini using the coupled channel Bethe-Salpeter (BS) 
equation. As shown in Refs. [2M31] using a certain renor- 
malization procedure, or in Refs. [35113^ using dispersion 
relations, in a certain partial wave one can write the BS 
equations in a factorized form, where the on shell potential 
and scattering amplitudes arc factorized out of the VGT 
integral term of the BS equations and, as a consequence, 
the integral coupled equations become simple algebraic 
equations. These BS equations are written as [26| 



T = [l-VG\-^V 



(1) 



where y is a 2 x 2 matrix accounting for the S-wave 
TTTT -^ TTTT, KK -^ KK and tttt — >■ KK potentials and 
G = diag(G'i, G2) is a diagonal matrix that accounts for 
the loop function of the two-meson propagators in the in- 
termediate states, given by 



|q|<<?n 



G, 



d^q 



(2^)3 2c^i(q)w2('Z) 



uJi{q)+bJ2{q) 



E'^ - [ioM) + W2{q)Y + ie' 



wi,2(q) = \ ■ml,2 + q 



(2) 



with meson masses mi , TO2 in channel j (in our particular 
case, / = 0, the masses mi and TO2 arc equal. To ease nota- 
tions, below we suppress the index j in the masses, etc.). 
We call KK channel 1 and tttt channel 2, and E = y/s 
is the center-of-mass (CM) energy. The function G re- 
quires regularization and we use, as in [26] . a cutoff regu- 
larization with a natural-size cutoff ^max = 904 MeV [55] , 
though it is also customary to use dimensional regulariza- 
tion with a subtraction constant to regularize the loops 
|36| . The choice of a particular regularization scheme does 
not, of course, affect our argumentation. Eq. ^ is just the 
relativistic generalization, consistent also with the meson 
statistics, of the integral of the ordinary non-relativistic 
Green's function 1/ {E — Hq) in the Lippmann Schwinger 
equation. By using the above cutoff, one obtains very good 
results for the tttt and KK scattering amplitudes up to 
about Vi - 1200 MeV [H]. 

We put now the same model in a finite cubic box of side 
length L and predict the discrete spectrum that emerges 
in such a box. Neglecting all partial waves except the S- 
wave [3, the only change to be made is to substitute the 
function G in Eq. ([TJ by G = diag (Gi, G2), where 



Gi — -TT 



|q|<9.r 



1 



UJi{q) + U2{q) 



i3 Z_. 2uji{q)u:2{q) E'^ - {ui{q) + uj2{q)f 
q = —n, n^lf , 

Lj 



(3) 



^ As rotational symmetry is broken in the box, partial waves 
mix in general, e.g. S-wave mixes with the L = 4 partial wave 
(G-wave). This effect is expected to be small and neglected in 
this study. 
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and again, uji^2(,q) = \/"^i2 + 9^- ^^ other words, the 
only change is to replace the integral over the continuous 
variable q in Eq. ([2]) by a sum over the discrete values, 
corresponding to periodic b.c. . 

The discrete spectrum of the system is given by the 
poles of the scattering T-matrix T{E) in a finite volume. 
It is interesting to show the relation of the present ap- 
proach to the one of Liischer for the case of one channel. 
To this end, assume first that the energy E is above thresh- 
old, and take into account the fact that the potential V 
is the same in a finite and in the infinite volume, up to 
exponentially suppressed terms. The poles emerge at the 
energies T~^{E) = 0, i.e.. 



2.2 Relation to the Liischer equation 

Next, we wish to demonstrate that Eqs. ([7]) and ^ are 
nothing but the ordinary Liischer equation in regions of 
energy and L where the Liischer approach is a good ap- 
proximation (see also Ref. [H]). To this end, let us use the 
identity 



1 



UJi +U2 



2 ui UJ2 E^ — (cji + ^2)^ + ie 

1 1 1 



1 



2E p2 
1 



1 



2uJi UJ2 ijJl + UJ2 + E 

1 1 



V-^{E)-G{E) =0. 



(4) 



4 wi W2 ti^i — a;2 — -B 4,101(^2 L02 — uJi — E 



(10) 



It can be easily seen that the above equation produces an 
infinite tower of discrete levels. 

Further, using Eq. (|3]), the scattering matrix in the 
infinite volume can be rewritten, for the discrete eigenen- 
ergies satisfying Eq. ([5]), as 



T{E) = {V-\E) G{E))-' = {G{E) G{E))- 



(5) 



Below, we shall use the normalization of Ref. [2S1 ■ In this 
normalization, the infinite- volume scattering matrix T{E) 
is related to the phase shift by 



T{E) = 



-8ttE 



pcot 6(j)) — ip 



(6) 



(below threshold, the ie prescription should be omitted). 
This identity should be used both in G{E) and G{E). 

Let us now take into account the fact that the singu- 
larities of the last three terms in the r.h.s. of Eq. pU]) are 
not located in the scattering region {E above threshold) 
and are separated from the threshold by a distance of the 
order of the particle masses mi^2- According to the regular 
summation theorem pTlf^ for such functions, the differ- 
ence between the sum and the integral is exponentially 
suppressed at large volumes. Consequently, the contribu- 
tion from the last three terms is exponentially suppressed 
in the difference G{E) — G{E) and can be neglected. This 
difference above threshold takes the form 



whevep = \^/'^{E'^, ml, ml)/{2E) a.ndX{x, y,z) standsioi G(E) — G(E) — < — V^ — / 
the Kiillen triangle function. Substituting this into Eq. (0) , \L^ ^-^ J 



(Pq 

(27r)3 



we arrive at 



pcoi5{p) = -^-kEI G{E) - G{E) 



_ip_\ 
8nEj 



above threshold 



(7) 



2Ep^ 



1 



. . l9l<9„« 

— — V 



2E L3 A^ „2 _ q2 
1 



. 2V \ *"ax + T, log — 



ip 

8ttE 



(11) 



Here, we suppressed the index "j" in the Green function, 
because the one-channel problem is considered. Note also 
that, for E above threshold, we get G{E) + ip/{8-KE) = 
ReG(i?), and the r.h.s. of Eq. ([7]) is real, as it should. 

Finally, we perform the analytic continuation of Eq. ^ 
below threshold. The l.h.s. of this equation is defined only 
above threshold. However, it is known that the l.h.s. obeys 
the effective-range expansion near threshold 



pcot S{p) = h - rp^ + Oip'^ 

a 2 



(8) 



where a and r denote the scattering length and the effec- 
tive range, respectively. Note also that this is an expan- 
sion in the variable p^ and hence the threshold is a regular 
point in this expansion. Performing now the analytic con- 
tinuation below threshold, where p = ij, we get 



where the ellipses stand for the exponentially suppressed 
terms. Below threshold, the last term in the r.h.s. gets 
replaced by —^/{8ttE). 

Moreover, as seen from Eq. (fTTjl . one may in fact re- 
move here the cutoff, sending (/max ^ 00. Indeed, one 
should obviously take a gmax such that p^ < g^^^ in the 
whole region of interest to us. If we sum and integrate from 
gmax to g^j^^, with q'^^^^ > g^ax, the denominator p^ „ q2 
is not singular and, according to the regular summation 
theorem, only exponentially suppressed corrections may 
arise. Finally, noting that (see, e.g., Ref. ^7\ ) 



11- 



l<?l<'?n 



^™ 1 T3 



E 



1 



p2 _ q2 



1 



1 1 2 
a 2 ' 



-8nE}^G{E) - [g{E) - ^ 
f below threshold 1 . 



27r3/2i 



27r2 
2oo(l,p2), 



P 



pL 

27r 



(12) 



(9) 



where Zqq stands for the Liischer zeta-function, we can 
directly verify that Eq. ([T]) (or its extrapolation below 
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threshold, Eq. ^) indeed coincides with the Liischer equa- 
tion 



c. X 27r 
p cot d (p) = — IT 
ij 



-3/2 



2oo(l,p') 



(13) 



and is cutofF-independent up to exponentially small cor- 

rectionfl 

Equation (|12[) can be used for practical purposes to 
evaluate the Liischer function, if one wishes, and as no- 
ticed in Ref. [37] its finiteness is explicitly shown. However, 
for practical purposes it is more convenient to keep the log 
terms in Eq. pip , which guarantees a faster convergence, 
and we obtain {A = gmaxi/(27r)) 



-2oo(l,P^) = Ihir 



1 



A 

E- 

n 



1 



F{p\A) 



F{p\A) = V47ryl 
[ (1 -|- 7rK/2 — K atan 1/k) , k — \J —p"^ j A else. 



VA — vl ^ for p ^ > 



(14) 



Below threshold, we have subtracted from G the analytic 
extrapolation of i Ini G. The convergence in practice is 
found for values of the dimensionless A as low as yl ~ p-|-3. 
The convergence is mostly limited by the artefacts intro- 
duced by the sharp cut-off A. It can be improved by av- 
eraging over several values oi A > p -I- 3. If one omits the 
log term in Eq. p4p like in Ref. [37], the same level of 
convergence is only reached for larger values of yl > j5-|- 6. 

We further comment on the question, whether the finite- 
volume corrections, stemming from the last three terms in 
Eq. (jlOp. are uniformly suppressed by an exponential fac- 
tor in the whole energy region measured on the lattice. 
From Eq. pH)) it is seen that, e.g., in the equal-mass case, 
the third and the fourth terms on the r.h.s. have a pole at 
-E = 0, so the convergence is not uniform in the vicinity of 
this point (note that the left-hand side of Eq. (jTU]) shows 
no singularity at E = 0). The singularity of the second 
term of the right-hand side is located further left in the 
energy plane. In general, for the case of non-equal masses, 
the singularities arise for the energies E < |mi — m2|. Of 
course, in case of one-channel scattering, this point lies 
below threshold, in the unphysical region. If the energy is 
taken larger, one may safely use the regular summation 
theorem and show that the contribution of the last three 
terms to the difference G — G is uniformly suppressed by 
an exponential factor. For smaller energies, this is no more 
true. Yet, for a case like ttZ', the singularity at |mi — TO2I is 
situated 2M^ below threshold and corrections due to the 
neglected terms show up much before the pole. One should 
also note that at energies above threshold, even if these ne- 
glected terms are exponentially suppressed, one still finds 



■^ Note that, despite the non-relativistic appearance of the de- 
nominator in Eqs. pip and p2|) that at the first glance looks 
like a finite-volume version of UCHPT with non-relativistic 
propagators [38II39| , the formalism of course possesses the rel- 
ativistic dispersion law. 



a contribution from those terms and their numerical rele- 
vance must be checked. This is done in appendix lAl 

In addition to this, one can readily see that there is 
a potential problem if one considers the multi-channel 
Liischer equation. Let m\ , mi and vn!-^ , m'l be the masses 
in the "heavy" and "light" channels, respectively, so that 
mi+m2 > m'l+mi- If |toi — TO2I > m'l+m'i, the singular- 
ity from the "heavy" channel comes into the physical re- 
gion of the hght channel and using the decomposition (10) 
can not be justified any more. Fortunately, it does not hap- 
pen in the case oi tttt — K K and nrj — KK coupled-channel 
scattering, which is considered in the present paper. 

Even if in the particular case of interest the redun- 
dant poles, which stem from the decomposition given in 
Eq. (10), do not sneak into the physical region, it is still 
legitimate to ask, whether the presence of such nearby sin- 
gularities may affect the numerical accuracy of the method 
(even above threshold, as we pointed out above). We would 
like to stress the following once more: the difference which 
we are trying to estimate, is exponentially suppressed at 
large L and thus beyond the accuracy claimed by the 
Liischer approach. However, as we can see in Fig. 1 be- 
low, the volumes for which the energies of the low-lying 
levels are of order of 1 GeV or less, are indeed not that 
large. 



2.3 Cutoff effects 



Finally, we wish to comment on the artifacts coming from 
the finite cutoff (/max- Keeping the sharp cutoff of natural 
size yields small unphysical discontinuities in the predicted 
energy levels. These discontinuities disappear if we use a 
smooth cutoff instead. 



G., 



Gs,: 



Gs.j + (Gj - Gs.j) 



dh 



FM) 



uji{q) + uj2{q) 



(2^)3 2uji{q) L02{q) E^ - (u^iq) + uj^iq))^ + le 



G 



s,j 



Here, 



-y- 

r 3 Z^ 



n\Q\) 



uji(q) +uj2{q) 



L3^2a;i(q)w2(g) E^ - (LUiiq) + LJ2{q)r 



F{q) 



(15) 



q" 



(16) 



We choose a — 1.2 GeV, m ~ 12. This form factor is 
almost equal to one up to approximately q ~ 800 MeV 
and then drops smoothly to zero at around g ~ 1.7 GeV. 
The only caveat in this procedure is that the form fac- 
tor F must be equal to one for the on-shell value of the 
momentum g, in order to obey unitarity exactly. In choos- 
ing F, we have demanded that it is different from one by 
less than l%c for all energies considered here. 
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Fig. 1. Energy levels as functions of the cubic box size L, 
derived from the chiral unitary approach of Ref. [2S] and using 
G from Eq. (US}. 
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Fig. 2. Solid line: The original tttt scattering S-wave phase 
shift 5q from the chiral unitary amplitude of Ref. [55] i Dashed 
line: the pseudophase extracted from the spectrum of Fig. (TJ 
using the one-channel equation for the level 2. 



3 Two-channel formalism 

3.1 General setting 

The generalization of the Liischer formalism to coupled- 
channel scattering is straightforward. The secular equa- 
tion that determines the spectrum is given by 

dct{t~VG) = l-VnGi-V22G2+{VnV22-V,\)GiG2 = 0. 

(17) 
The energy levels which have been determined from this 
equation, using the input potentials Vij from Ref. |26j , are 
shown in Fig. [T] 

In Refs. PHHS], we have introduced a very useful quan- 
tity called pseudophase. We shall often refer to this quan- 
tity below, because it is convenient in the discussion of the 
Liischer equation. By definition, the pseudophase is the 
phase extracted from a given energy level, using the one- 
channel Liischer equation ([7]) (i.e., neglecting the channel 
with a higher threshold). From the definition, it is clear 
that, well below the inelastic threshold, the pseudophase 
coincides with the usually defined scattering phase, up 
to the exponentially suppressed finite- volume corrections. 
Above the inelastic threshold, a tower of resonances emer- 
ges in the pseudophase, which reflect the opening of the 
second channel. 

In Fig. [21 we show the tttt pseudophase, extracted from 
the energy levels in Fig.[TJ It is seen that, below 900 MeV, 
the pseudophase is very close to the tttt scattering phase 
Sq. Above 900 MeV, the pseudophase rises faster than 5q, 
owing to the presence of the KK threshold. 



3.2 Multi-channel approach 

Next, we turn to the central problem for the two-channel 
(generally, the multi-channel) approach. As mentioned al- 
ready, unlike the one-channel case, a single measurement 



of a value of the energy E{L) at a given L can not provide 
the full information about three independent quantities 
Vij{E{L)) at the same energy. The variation of L does 
not help, because the energy also changes. In order to 
determine all Vij independently (and, thus, the pole posi- 
tions, which are eventually determined by the Vij), there 
are several options: 

i) In Ref. [5S] we have proposed to use twisted b.c. for 
the s-quark. If with Q < 6i < 27r, z = 1, 2, 3 denotes 
the twisting angle, imposing the twisted b.c. is equiva- 
lent to adding ztO/L to the relative three-momentum 
of a kaon (anti-kaon) in the CM frame, whereas the 
relative momenta in the tttt intermediate state do not 
change (Below, for simplicity, we shall always use the 
symmetric choice 9i = 6). Consequently, the expres- 
sion of the TTTT loop function does not change, whereas 
the KK loop function becomes 6'-dependent (the per- 
tinent expressions can be found in Ref. [25] V Such a 
procedure is very convenient, because it allows to move 
the KK threshold and thus perform a detailed scan of 
the region around 1 GeV, where the /o(980) is located. 
The idea behind using twisted b.c. is that in this case 
one has two "external" parameters L and 9, which can 
be varied without changing the dynamics of the sys- 
tem. In particular, one may adjust these so that a given 
energy level does not move. On the other hand, since 
the loop functions, entering the secular equation (fT7|) . 
depend on L and 9, for different values of these param- 
eters one gets a (non-degenerate) system of equations 
that allow the determination of the individual matrix 
elements Vij at a given energy. 

For illustrative purposes, in Fig. [3] we display the 9- 
dependent spectrum for 9 = (periodic b.c), 9 = 
tt/2 and 9 = n (antiperiodic b.c), obtained by using 
UCHPT in a finite volume. It is seen that the spec- 
trum is strongly sensitive to twisting at the energies 
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Fig. 3. Solid lines: Spectrum E{L) as in Fig.[T] Dashed lines: 
Applying twisted b.c. = n as in Ref. [SS], only for the KK 
channel and for the levels 2,3, and 4; Dashed dotted lines: with 
e = 7r/2. 



around 1 GeV, which is exactly the energy region we 
are interested in. 

We realize that it could be quite challenging to im- 
plement this idea (with a complete twisting, including 
the sea quarks) in present-day lattice simulations. A 
partial twisting (of only the valence quarks) can be 
performed more easily. The question, whether such a 
partial twisting will still enable one to extract useful 
information about the properties of the multi-channel 
resonances, is still open. We plan to address this (very 
interesting) question within the EFT framework in the 
future, 
ii) Introducing the twisting angle 6 gives us an additional 
adjustable parameter that allows one to determine the 
individual matrix elements Vij . The same goal can be 
achieved without twisting, performing the calculations 
on the asymmetric lattices L^ x Ly x L^ where, gen- 
erally, Lx^y^z are not all equal Q Note that scattering 
in an asymmetric box has already been studied on the 
lattice, see Refs. [1011111111113. 

Suppose Lx.y^z have been adjusted so that the value 
of the energy E{Lx^y^z) = E stays put. In practice 
one evaluates the trajectories of the levels, and cuts 
them by the line of £' = const to determine the Li. 
The loop functions Gj depend both on the energy E 
and the box configuration given by L^^y^z- Note that 
L^ in Eq. psp is now replaced by L^LyLz and qi = 
2Tmi/Li, i = 1,2,3, rii £ Z. We label the different 
configurations with the same energy by an index a = 
1, 2, • • • , and denote the pertinent values of the loop 
functions by G^""^ (E). Performing the measurement for 
three different configurations and solving the secular 
equation, one can determine the matrix elements of 
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the potential separately, 



V,%-VnV2: 



Gr 



(18) 



where 



qW nW ^(i)/=.(i)' 



G 



Gl = GY' G'^'> gY'g 



^(2) 



1 2 

i(2)/=,(2) 



(19) 



i(3) M3) 



{3)M3) 



Uj U2 Uj U2 



Effects of additional partial wave mixing are neglected here. 



In Fig. IHwe show a particular energy level in the vicin- 
ity of 1 GeV, calculated for different configurations of 
the asymmetric box. It is seen that the level is not 
very sensitive to the configuration. For this reason, a 
larger error bar is a priori expected, if observables are 
extracted from these data. 

For illustrative purposes, in Fig. [51 we display different 
physical observables in the tttt and KK scattering, that 
can be extracted from the energy levels by using the 
above-mentioned methods, no errors yet assigned. As 
seen, both methods have in principle the capability to 
address the problem. 

For completeness, one should mention that measuring 
at least 3 different excited levels at the same energy but 
different values of L could, in principle, also achieve the 
goal. Such a proposal, however, looks quite unrealistic in 
the view of the present state of art in the scalar meson 
sector. 



3.3 Threshold effects 

At this place, we would like to comment on the threshold 
effects, which have a capability to seriously complicate the 
search of the near-threshold resonances - exactly the task 
we are after. The resonances on the lattice reveal them- 
selves in the form of a peculiar behavior of the volume- 
dependent energy levels near the resonance energy called 
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"avoided level crossing" [?3|. For example, in Fig. [31 the 
solid curves (levels, corresponding to the periodic b.c), 
display a perfect avoided level crossing slightly below 1 
GeV. Does this signal the presence of the /o(980)? 

In order to answer this question, it suffices to change 
the parameters of the potential slightly, making the /o(98G) 
disappear. This can be achieved, for example, by merely 
reducing the strength of the potential Vn that describes 
the scattering in the KK -^ KK channel. Namely, re- 
placing Vii — > ryyii and varying rj between 1 and 0, it is 
seen that the /o(980) disappears already at ?7 ~ 0.8 (the 
pole moves to a hidden sheet). Below, we shall study the 
dependence of the energy levels on the parameter rj. 

Let us first stick to the periodic b.c. . In the left panel 
of Fig.[Slone immediately sees that the avoided level cross- 
ing in the energy levels persists even for vanishing 77, when 
there is no trace of the /o(980) any more in the scatter- 
ing phase (see Fig. [TJ. This shows very clearly that the 
avoided level crossing in this case has nothing to do with 
the /o(980) and is due to the presence of the KK thresh- 
old. The modified problem has memory of the KK chan- 
nel because the V12 term is still kept and connects the two 
channels. 

The pseudophase which is displayed in the left panel 
of Fig. \7\ also demonstrates this statement. We see that 
it barely changes with 77 — > 0, whereas the true scatter- 
ing phase is subject to dramatic changes. A rapid change 
of the pseudophase by tt around 1 GeV is related to the 
presence of the KK threshold and not to the /o(980). 

Now, the essence of the problem becomes clear. It 
seems that the presence of the threshold masks the pres- 
ence of the resonance, and the energy levels become largely 
insensitive to the latter. How accurate the lattice data 
should be in order to provide a reliable extraction of the 
resonance parameters in the vicinity of a threshold? 

It should be pointed out that the use of twisted b.c. 
gives one possible solution to the above problem. Namely, 
as already mentioned above, as the twisting angle 9 in- 
creases, the KK threshold moves up, whereas the /o(980) 
stays put. Consequently, for larger values of 6, these two 
are well separated, and the pseudophase reveals a larger 
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Fig. 6. Energy levels, different values of 77. Left panel corre- 
sponds to periodic b.c. 9 = 0, right panel to twisted b.c. 9 = n. 
Solid lines: 77 = 1; Dashed lines: 77 = 0.8, Dash-dotted lines: 
ri = 0. 
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Fig. 7. Thick lines: true tttt scattering phases, different values 
of rj; Thin lines: pseudophases, different values of 77. Left panel 
corresponds to periodic b.c. 9 = 0, right panel to twisted b.c. 
= n. Solid lines: 77 = 1; Dashed lines: 7; = 0.8, Dash-dotted 
lines: 7; = 0. 
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dependence on the parameter ?y. This is demonstrated in 
the right panels of Figs. [S] and [71 

It remains to be seen, whether the use of twisted b.c. 
could be advantageous for resolving near-threshold reso- 
nances. Below, we address this issue by means of fits to 
synthetic lattice data. 



3.4 The framework based on UCHPT 

Finally, we would like to turn to the framework of UCHPT 
in a finite volume, which is a central issue considered in the 
present paper. In section \J72\ we did not assume any spe- 
cific ansatz for the potential Vij which is being extracted 
from the data. It is clear, however, that, changing the box 
geometry and/or the twisting angle, we can not arrive 
exactly to the same energy, and some kind of an interpo- 
lation procedure will be needed. This could be achieved 
by evaluating many points belonging to the same level for 
different values of L and performing a best fit to the data 
which will return an accurate trajectory with smaller er- 
rors than the individual ones (see, e.g., Refs. [iSlHSlITT] ) Fl. 

Alternatively, one may choose the parameterization of 
the potential Vij{E) as in UCHPT (e.g., as in Eq. 1^ 
below) , and then fitting the parameters of Vij to the data 
on the lattice introduces certain model-dependence, but 
provides a fair input based on the successful UCHPT that 
allows the interpolation of the lattice data. With more 
lattice data available, Vij as functions of the energy can 
be fitted more accurately, making the initial ansatz for 
it less relevant. We would like to stress that the present 
ansatz for the potential, in our opinion, is simple, rather 
general and theoretically solid. As mentioned above, it has 
been successfully used in the past to accurately describe 
the scattering in the mr/Trr] and KK systems (at physical 
values of the quark masses), see Refs. pSl[?7l[^ [^. This 
means that the effect of the higher-order terms in the po- 
tential should be small. Of course, this assumption could 
be checked by incorporating such higher order terms as 
it has been done e.g. in the case of K~p scattering [55]. 
It can be expected that, in the first approximation, the 
functional form of the potential stays intact when we go 
to higher quark masses that are used in present lattice 
simulations (the numerical values of the parameters may 
depend on the quark masses). 

Below, performing the fits to the synthetic lattice data 
with errors, we shall explicitly demonstrate how the ap- 
proach works. 



4 Fit to the energy levels 

4.1 Twisted boundary conditions versus asymmetric 
boxes 

Below, we shall carry out a preliminary numerical esti- 
mate of the efficiency of the two different schemes which 
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were discussed above. To this end, first we produce the en- 
ergy levels by using UCHPT in a finite volume. Then we 
assume that each data point E{L) is produced with some 
uncertainty AE{L). Reversing the dependence on L, one 
obtains L{E) and AL{E) (in case of an asymmetric box we 
have taken Li ~ ai L for three different choices of ai and 
we take the same AL as in the symmetric case) . In general, 
for a given E{L) and AE{L), the inverse AL{E) is a very 
complicated function of energy, depending on the tangent 
of the curve. Since we are making only a rough estimate 
here, we simplify our task by assuming a constant error 
AL = 0.02A/^^. Further, using the von Neumann rejec- 
tion method, we generate a Gaussian distribution for the 
variable L centered around the exact solution a = L{E) 



f(L) 



1 



ALV2^ 



exp 



JL-a)' 
2{ALf 



(20) 



* We are aware that present day lattice simulations do not 
provide more than a few volumes. 



We perform the calculations of the infinite-volume phase 
shifts for each value of L and determine the average value 
and uncertainty, using Eq. (fT5|l . The results for two dif- 
ferent methods are shown in Fig. [S] Below and above the 
KK threshold, both methods deliver very similar phase 
shifts with a comparable error. Around the KK thresh- 
old, as expected, the use of twisted b.c. produces phase 
shifts with an error smaller than in the other method. 



4.2 Extraction of the /o(980) pole position from the 
data 

After the preliminary study of the different schemes car- 
ried out in the previous subsection, we turn to the main 
topic of our paper and perform the fit of the data by using 
UCHPT in a finite volume. To this end, we consider the 
procedure for producing various synthetic data sets and 
the strategy of the fit. 

The central values £'„(L), n = 0, 1, • • • are produced 
by solving the secular equation \YJ\ with the potentials 
from UCHPT, by using periodic/antiperiodic b.c. (below, 
we do not study the case of the asymmetric boxes) . At the 
next step, a constant error AE =10 MeV is assigned to 
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each of the data points. We consider three different data 
sets in the vicinity of the /o(980), each of them containing 
13 data points (see Fig:[n|): 

1. The data set from the levels 2 and 3, using periodic 
b.c. 

2. The level 2, using both periodic and antiperiodic b.c. 
As we shall see, the use of the twisted b.c. enables one 
to achieve a better accuracy than in set 1. 

3. The levels 2 and 3, periodic b.c, modified input (the 
parameter 77 = 0.6, so the /o(980) pole is absent). As 
seen from Fig. [HI there is no qualitative change of be- 
havior of the energy levels in sets 1 and 3. It is im- 
portant to check, whether the method is capable to 
distinguish between these two cases. 

Further, for each set the data Ei{L) shown in Fig. [5] 
are fitted by using linear functions in the potential Vij , 



V,. 



Mj+hj{s-AMl) 



(21) 



which is used to calculate E[. 



(theory) 



(i) by using Eq. P7)) . 
The six parameters of the potential a^j and bij are as- 
sumed to be free and are determined from the fit to the 
spectrum. Note that, since the input potential which was 
used to produce the data points is also linear in s, the 
X^-value of the best fit is equal to 0. 

To obtain the error on the extracted quantities (phase 
shifts and pole positions), random Cy, bij are generated 
within the limits given by the parameter errors, that are 
determined by the previous fit. Of these combinations, 
only those are accepted for which the resulting y^ is smaller 
than Xbest + 1- F'o'^ each such event, the tttt phase is evalu- 
ated and the pole position in the complex energy plane is 
determined. The resulting bands for the phase shifts and 
extracted pole positions are shown in Figs. [TUl and [TT] 

From these figures one observes that, for all fits to 
the synthetic data with /o(980) (sets 1,2), which obey the 
Xbest "I" 1 criterion, the phase shift shows a resonant behav- 
ior and a pole is found in the complex plane on the first 
KK sheet (second sheet for the tttt channel). Also, in all 
fits to data without /o(980) (set 3), the phase shifts show 
no resonant behavior, and no pole is found in the vicinity 
of the KK threshold — the poles are in this case on a 
hidden sheet far from the physical region. Thus, with the 
method proposed here, it is indeed possible to distinguish 
clearly the presence or absence of a resonance at the KK 
threshold. 

Moreover, as the figures [TO] and [11] show, using 13 
data points with a 10 MeV error, the pole position of 
the /o(980) can be determined quite precisely. As seen in 
Fig. [TT] the data set 2 gives a more precise determination 
of the pole position than the data set 1: it is advantageous 
to use data from a lower level but with different bound- 
ary conditions, than to use the data from higher excited 
levels. 

If fewer lattice data are included in the analysis than 
those of Fig. [5] the six-parameter fit is less constrained. 
Then, the spread in extracted pole positions and phase 
shifts can increase. This is reflected in a blow-up of the 
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10 



M. Doring et al.: Unitarized Chiral Perturbation Theory in a finite volume: scalar meson sector 




850 



900 



950 
E [MeV] 



1000 1050 1100 



Fig. 10. Extracted phase shifts corresponding to the three sets 
shown in Fig. [S] Dashed lines: The calculated phase shifts by 
using the approach of Ref. [^Hi; Filled (cyan) bands: Recon- 
struction of these phase shifts. 



bands shown in Fig. [HI in regions without data. Then, 
adding a data point precisely in this region helps constrain 
the fit. In an actual lattice calculation, such observations 
could help to guide the selection of lattice sizes L for which 
the levels are calculated. 



4.3 Statistical scattering of the data (Set 2') 
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Fig. 11. Extracted pole positions. Filled triangle: original pole 
position [2S]; Shapes: Pole positions reconstructed from the 
data sets 1, 2, and 2', see Fig. [§1 For the data set 3, the poles 
are found either on the hidden sheets, or very far from the KK 
threshold. 



distribution with (7 = 5 MeV for the displacement of the 
centroids. Using Eq. (PU|) for every data point, 13 new 
points for Set 2 are generated and then fitted. The result- 
ing fit leads to a slightly modified pole position and phase 
shift. To obtain, however, reliable results, this procedure 
has to be repeated multiple times. The multiple generated 
data sets are called Set 2' in the following. To obtain con- 
fidence regions for levels and phase shifts, the multiple fits 
to Set 2' are analyzed. For example, for the resulting phase 
shifts, at every energy E the mean value and standard de- 
viation for all fits are calculated. The resulting band is 
indicated with the dashed lines in Fig. [TUl The same can 
be made for the levels themselves as indicated with the 
dashes lines in Fig. [21 

To estimate the combined effect of uncertainty from 
the previously discussed Xbest + 1-criterion and displaced 
centroids, we calculate the average x^^^, from the multiple 
fits to set 2', with the result x^ve ~ !•'''■ The overall uncer- 
tainty is then estimated from applying a X^est + Xavc + 1" 
test on set 2, i.e. the data set without displaced centroids. 
The resulting bands are indicated with the hatched areas 
in Figs. [Hi and [TUl As expected, their widths is approx- 
imately given by the sum of widths from the Xbcst + 1" 
criterion [solid (cyan) areas] plus the 1-cr region of the 
multiple fits to set 2' [dashed lines] . In Fig. [TTJ the com- 
bined uncertainty for the pole position is shown [hatched 
area] . In summary, we observe that the scattering of the 
centroids leads to a moderate increase of the uncertain- 
ties for levels, phases and pole positions. With a Gaussian 
distribution of 5 MeV width for the displacement of the 
centroids, the uncertainties induced by 10 MeV error bars 
are increased by maximally 50 %. 



In an actual lattice simulation the finite statistics will not 
only contribute to the error bar of a data point but also 
lead to a displacement of its centroid. To study this ef- 
fects, we concentrate on data set 2 and choose a Gaussian 



4.4 Systematic uncertainties 

A comprehensive treatment of systematic uncertainties 
would require to address the actual lattice action that 
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energy spectrum (right panel). The solid curve and the circles 
correspond to dj = 0. The dashed line on the left-hand side 
and the squares on the right-hand side correspond to c^ 7^ 
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generates the levels, but that is far beyond of what can be 
possibly done in the present framework. What can be stud- 
ied, however, is the scenario when the chosen fit potential 
cannot achieve a perfect fit of the lattice data. As men- 
tioned above, the choice of the potential of Eq. (PT|) repli- 
cates the lowest-order chiral interaction that has been used 
to generate the synthetic data. To test the more realistic 
situation, in which the lattice data can not be perfectly 
fitted, apart from the statistical scattering of the data dis- 
cussed in the previous section, we introduce a small term 
~ Cij (s — 4M|-)^ in the potential that is used to gener- 
ate the synthetic data. The constants c^ are chosen to be 
cii = C22 = — c, C12 = C21 = +c where c = 310.8 GeV~ . 
The effect of introducing this term on the phase 5^ and on 
the energy spectrum (periodic b.c.) is seen in Fig.[T21 The 
fit to the modified data points (squares in Fig. [T^ and 
subsequent analysis are then performed using the same 
potential as before, i.e. Eq. (pijl . 

The resulting x^ost does not vanish any more, and 
the extracted pole positions and phases are systemati- 
cally shifted as shown in Fig. [121 case (a) and (b) (the 
latter correspond to data sets 1 and 2 in Fig. [5] In addi- 
tion, to make these areas distinguishable, we have assigned 
an error AE = 5 MeV in case (a), instead of 10 MeV.). 
This kind of systematic error is inherently tied to the as- 
sumption made on the functional form of the potential 
(cf. Eq. (HI])): one could, of course, introduce higher or- 
der terms in the potential to extract phase shifts and pole 
positions. Then, Xbost would be close to zero again, but 
the fit would be much less constrained given the higher 
number of free parameters. The spread of pole positions 
and extracted phases would then immediately increase as 
compared to the present results. Still, it would be worth 
to quantify the effects of the higher order terms on the 
accuracy of these determinations in future studies. 

Apart from the discussed higher order terms for the 
potential V, another source of uncertainty is given by the 
cutoff (7max that is used to determine phase shifts and pole 
positions once the Vy are fitted: the term G of Eq. p7p 
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depends on this (jmax- To test the dependence, we have 
extracted the pole positions with gmax = 755 MeV instead 
of using the standard value of 904 MeV that has been used 
to generate the synthetic data. The result is shown as case 
(c) in Fig. [121 Obviously, the functional form of Eq. (PT|) in 
Vij is fiexible enough to compensate for this large change 
in gmax and as a result the pole position is shifted by just 
a few MeV. Systematic effects from the cut-off dependence 
arc, thus, very small. 

We have made another test which can be of use if one 
could make a rough estimate of the size of the systematic 
uncertainties of the lattice results. Systematic uncertain- 
ties would most likely produce deviations of the energy 
levels in the same direction. We have assumed a constant 
shift of the energy levels, increasing them by 5 MeV, and 
we have found a decrease of about 5% in the phase shifts, 
and viceversa. 

To summarize, a careful examination of all kind of er- 
rors has been carried out in this subsection, illustrating 
the limitations of the proposed method. Given only a few 
lattice data with large error, UCHPT in finite volume is 
certainly a useful tool to make reliable statements on the 
presence or absence of resonances, and even to give a good 
first estimate on phases and pole positions. 



4.5 The ao(980) 

The situation concerning the ao(980) in UCHPT is in- 
triguing. For example, in Ref. [55], the ao(980) is gener- 
ated from the lowest order chiral Lagrangian in the irrj and 
KK channels. At present, there exists no consensus about 
its nature and exact location. For example, changing the 
cutoff (jmax in the approach of Ref. [SHj, one may easily 
move the pole below/above the KK threshold. Namely, 
for_gniax = 1094 MeV, the pole of the oq is below the 
KK threshold and one observes the full resonance shape 
(see the left panel of Fig. [H)) . For lower cut-offs, say 
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all. Note also that, as_seen from Fig.[T51 the small efFeet of 
the sub-threshold KK channel is visible in the extracted 
phase for the highest energies above 600 MeV. Also, for 
the same reason, there is an off-set between the true pole 
position (dark circle, right-hand panel) and the extracted 
ones (large shaded area). To verify this, we have generated 
pseudo data using only the TTTr-channel in the hadronic 
model |26] . Then, extracting pole positions with the one- 
channel formalism, the resulting ellipse should have the 
true pole position (of the TTTr-reduced hadronic model) in 
its center. This is indeed the case as has been checked. 

The /o(600) (or a) has been predicted in most UCHPT 
approaches (see, e.g., Refs. pSll^HU] ). Recently, its exis- 
tence has been rigorously proved by using ChPT com- 
bined with Roy equations |54| . It would be intriguing, if 
the already existing lattice studies of the lowest scalar res- 
onance in QCD p fTTllTB] could be supplemented by the 
novel method of the analysis based on UCHPT, in order 
to facilitate the extraction of the pole position in the com- 
plex plane. 



Qmax — 904 MeV and g^ax = 630 MeV, the pole disap- 
pears from the second Riemann sheet and instead of a 
resonance shape, a pronounced threshold cusp structure 
is visible. This was already observed in Refs. P^49j . 

This is an interesting situation, because there are sev- 
eral examples in the literature where it is not clear if an 
observed structure is a resonance or rather a cusp (exam- 
ple: a disputed pentaquark in 772 -^ r\N [5U1[ST1I5^1[^ ) . 
Using UCHPT in a finite volume, one may predict the 
dependence of the energy levels on the cutoff parameters 
and see, whether there is a clear-cut distinction between 
the situations with a resonance or a cusp. 

In order to show how the different scenarios influence 
the energy levels in the box, we carry out the calculation 
and show the results in Fig. [TH In the left panel of this 
figure we show the lineshapes in -kt] — > -kt] for three differ- 
ent cutoffs. In the right panel, we display the energy level 
1 for the same values of the cutoff, both for periodic and 
antiperiodic b.c. As one observes from this plot, the level 
is quite sensitive to the value of the cutoff. Consequently, 
it is expected that the extraction of the pole position from 
future lattice data can be performed accurately, resulting 
in a clear-cut resolution of the resonance/cusp scenario. 



5 Conclusions 

In this paper, we discuss the extraction of the parame- 
ters of the /o(600), /o(980) and ao(980) resonances from 
lattice data. In order to facilitate this extraction, we use 
UCHPT in a finite volume in the fit. Fitting synthetic 
data, we have demonstrated that the approach works, and 
the pole position for the above resonances can indeed be 
extracted by analyzing the volume-dependent energy spec- 
trum in the vicinity of the resonance energy, provided suffi- 
ciently many volumes are simulated. The different sources 
of errors are analyzed in detail. The key point is that the 
use of the phenomenological input from UCHPT stabi- 
lizes the fit. It is, thus, very challenging to apply these 
results to present and forthcoming lattice data in order 
to extract valuable information concerning hadronic reso- 
nances. From the theoretical point of view one should also 
keep in mind that in field theory there are parts of the 
kernel of the Bethe-Salpeter equation which are volume 
dependent, although exponentially suppressed. Quantita- 
tive studies of these terms and the extent of the large L 
suppression would be a very good complement to the work 
we have carried out here. 



4.6 The /o(600) 

It comes rather as a surprise that the procedure seems to 
work even in case of the /o(600). In order to do this, we 
have resorted to the 1-channel Liischer formalism, because 
the energy is well below the KK threshold. Only the level 
1 with periodic b.c. was fitted (see Fig. [13 left panel). The 
fit resulted in the phase shifts shown in the center panel 
and in the pole positions displayed in the right panel of 
the same figure. As seen from the figure, the spread in 
pole positions is larger than in case of /o(980), but, in our 
opinion, it is remarkable that one is able to address the 
question of the extraction of the /o(600) on the lattice at 
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A Quantitative comparison to the Luscher 
approach 

As we have seen in Sec. 12.21 the connection of our approach 
with the Luscher approach conies through the replacement 



GjiE) -^Gja{E), where 
G,a{E)-G,{E) = 



in^/^EL 



2^00 (1;f 



^2\ 



WttE^ 



\^/\E\mlml). (22) 



It is straightforward to see that, substituting G in Eq. ([T|) 
by GjA of Eq. (P^ . one arrives at the two-channel Luscher 
equation from Ref. [25| . which is written in terms of the 
A'-matrix. 

Next we discuss the differences of the functions Gj (E) 
of Eq. P^ and Gja (E) . In Fig. [TB] we show the results 
for two different values of L for the tttt and the KK 
channels. We can see that for the tttt channel the differ- 
ences for L = 1.7 Af^^ are rather large for energies below 
250 MeV, which is already below the tttt threshold. Even 
above threshold, the differences are sizable. These differ- 
ences are smaller for a bigger L, L ~ 3.0 M~^, as expected 
for exponentially suppressed terms. In both cases the func- 
tion Gja{E) has a pole at -E = 0, while Gj{E) is finite 
there. For the case of the KK channel, the trend is simi- 
lar to the TTTT channel, but the region where Gja{E) and 
Gj (E) are practically equal is larger than in the tttt case. 

It is interesting to see, how big the effect of this dif- 
ference is on the energy levels in the box. This is seen in 
Fig. [13 left panel. As we see from Fig.[T71 left panel, if L > 
1.6 M~^, the corrections to the energy levels 2, 3, ...are 
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Gja{E) of Eq. H22|) [dashed lines], for the tttt and KK channels. 
Two cases for L — 1.7/Af,r and L = 3.0/A/^ are shown. 



completely negligible. Therefore, we concentrate on level 
1 where the difference between the solid and dashed curves 
can not be neglected. We calculate the tttt pseudophases, 
see Fig. [T71 central panel, and the position of the poles of 
the /o(600) corresponding to these two curves, see Fig.lTTl 
right panel 0. 

As one sees from Fig. [T71 the pseudophases are fairly 
the same until 500 MeV, from where they start diverging. 
At L = 2M-1 that corresponds to ^ ~ 700 MeV, the 
difference between two pseudophases is around 5 degrees. 
If the energy grows (the volume decreases), the difference 
between the two pseudophases becomes, as expected, big- 
ger. 

In Fig. [T71 right panel, we display the shift of the pole 
positions due to the effect of using the relativistic prop- 
agator. The same method as in section 14.61 was used to 
produce the ellipses. Case (a) shows the pole positions ex- 
tracted from levels generated using the relativistic propa- 



For small volumes, formally exponentially small terms 
exp(— n Mtt L) are no longer suppressed and would have to be 
evaluated, too. 
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gator (identical case to the one shown in Fig. [T5)) . In case 
(b) the pole positions have been extracted from the level 
generated with Gja{E) of Eq. (H2), i.e. from the dashed 
curve shown in the left panel. As we see, the real part of 
the pole undergoes a shift of order of 10 MeV, whereas 
the width (twice the imaginary part of the pole position) 
changes by approx. 40 MeV. 

As we see, the uncertainty of the method which stems 
from the error of 10 MeV in energy attached to each data 
point, is still larger than the change induced by the use of 
the relativistic propagators. Should one aim at an accu- 
racy of better than 40 MeV in the width, first, the errors in 
the lattice data should improve beyond 10 MeV and, sec- 
ond, from our study it follows that the standard Liischcr 
approach is not accurate enough in this case. 

So far, we were concerned how a shift in the level trans- 
lates into a shift of the pseudophase and pole position; one 
may also ask what happens if Gja{E) of Eq. (22) is used 
to extract the pole position from a given level, instead of 
using the relativistic propagator Gj (E) of Eq. ([TS]) as done 15. 
before. The effect is shown as case (c) in the right-hand 
panel of Fig. [T71 Compared to case (a), one observes a 16. 
systematic shift of the results, which is of similar size as 
the previously discussed case (b). 
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